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A two-phonon version of the relativistic quasiparticle time blocking approximation 
(RQTBA-2) represents a new class of many-body models for nuclear structure cal- 
culations based on the covariant energy density functional. As a fully consistent ex- 
tension of the relativistic quasiparticle random phase approximation (RQRPA) , the 
two-phonon RQTBA implies a fragmentation of nuclear states over two-quasiparticle 
and two-phonon configurations. This leads, in particular, to a splitting-out of the 
lowest 1~ state as a member of the two-phonon [2 + £g>3 _ ] quintuplet from the RQRPA 
pygmy dipole mode, thus establishing a physical mixing between these three modes. 
The inclusion of the two-phonon configurations in the model space allows to describe 
the positions and the reduced transition probabilities of the lowest 1~ states in iso- 
topes 116 ' 120 Sn as well as the low-energy fraction of the dipole strength without any 
adjustment procedures. The model is also applied to the low-lying dipole strength 
in neutron-rich 68 > 70 > 72 Ni isotopes. Recent experimental data for 68 Ni are reproduced 
fairly well. 

The theoretical description of nuclear low-lying dipole strength remains among the most 
important problems in nuclear structure and nuclear astrophysics [l| • Measurements of the 
dipole strength by means of high resolution nuclear resonance fluorescence 
resolve the fine structure of the spectra below the neutron threshold. Unique spectroscopic 
information about neutron-rich medium-mass and heavy nuclei have been obtained in recent 
experiments with Coulomb dissociation j^, 9] and virtual photon scattering [10]. This offers 
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exciting opportunities for microscopic nuclear structure models to describe the fine structure 
of the dipole spectra below the neutron threshold and to help in analyzing the experimental 
data. 

Measurements of the low-lying dipole strength above and below the neutron threshold 
are performed with different nuclear reactions, which have reduced sensitivity in the area 
around the threshold. Therefore, a correct comparison of the calculated pygmy strength 
with the data is still problematic. The collectivity of the pygmy mode is another subject of 
discussions. Self-consistent relativistic QRPA calculations produce a hi ghly collective 
pygmy mode, in contrast to the results of the non-relativistic approaches 12j, ll3|. However, 
in Ref. 14| it is pointed out that the collectivity of the pygmy mode is restored within 
non-relativistic QRPA calculations with Skyrme forces when a fully self-consistent scheme 
is employed. 

Since the pygmy dipole mode has essentially surface nature, it mixes with other surface 
modes, especially with low-lying ones. This supposition has been confirmed by explicit 
RQTBA calculations 15, Ji| for tin and nickel isotopes and N=50 isotones. It has been 
found that the pygmy mode, arising in RQRPA as a single state or as very few low lying 
dipole states with isoscalar character, is strongly fragmented over many states in a broad 
energy region due to the coupling to phonons. As a result, some fraction of the strength is 
located well below the original position of the RQRPA pygmy mode. Compared to existing 
data, the RQTBA describes the total fraction of the dipole strength below the neutron 
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threshold very reasonably [15J, H6( . However, in order to account for the fine structure of the 
spectrum, more correlations should be included in the microscopic model. In particular, it 
has to reproduce the lowest dipole state in vibrational nuclei, which is identified as a member 



of the quintuplet 2f ® 3j_, as it was predicted in Refs. [171, [l8| and observed in spectra of 



spherical nuclei 
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In the present work, the two-phonon version of the quasiparticle time blocking approxi- 
mation (QTBA) proposed in Ref. 23( is employed to introduce correlations between the two 
quasiparticles within the 2q®phonon configurations of the RQTBA. Therefore, a fragmen- 
tation of the nuclear states over the two-phonon configurations appears in the excitation 
spectra in addition to the spreading over the two-quasiparticle states. As well as in the con- 
ventional RQTBA, excitations in Fermi-systems with even particle number are described by 
the Bethe-Salpeter equation (BSE) containing both static and energy- dependent residual in- 
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FIG. 1: The correspondence between the 2q(g)phonon amplitude of the conventional phonon 
coupling model and the two-phonon amplitude of the two-phonon model in a diagrammatic 
representation. Solid lines with arrows and latin indices denote the single-quasiparticle nucle- 
onic propagators, wavy curves with greek indices the phonon propagators, empty circles represent 
phonon vertices, and grey circles together with the two nucleonic lines denote the RQRPA transi- 
tion densities (see text). 

teractions. This equation for the nuclear response function R(u) in the doubled quasiparticle 
space reads: 

kf,k(,ri" 

where the indices ki run over single-particle quantum numbers including states in the Dirac 
sea and the indices rj,rj',T]" numerate forward (+) and backward (— ) components in the 
doubled quasiparticle space. The quantity R^Z ( u ) ~ V ~~ ~~ Ek 2 ) describes the free 
propagation of two quasiparticles with their Bogoliubov's energies E kl and E^ 2 in the mean 
field between subsequent interactions with the amplitude 

w&Mkfr) = VStw* + K*,jw>) - Kk^M^'- ( 2 ) 

Here V is the static par t of the effective residual quasiparticle interaction. In the present 



work, as in Refs. 



15 



16|, it is derived from the covariant energy density functional (CEDF) 



with the parameter set NL3 24| as a one-meson exchange interaction with a non-linear 
self-coupling between the mesons. Pairing correlations are introduced into the relativistic 
energy functional as an independently parameterized term. In the 2q®phonon version of 
the RQTBA the energy- dependent residual interaction <&(u) is derived from the energy- 
dependent self-energy by the consistency condition and calculated within the quasiparticle 
time blocking approximation. This approximation means that, due to the time projection 
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in the integral part of the BSE, the two-body propagation through states with a more 



23| 



complicated structure than 2qCS>phonon is blocked 

In the present work we go a step further and consider a modified version of the RQTBA 
which includes additional correlations between two quasiparticles inside the 2q®phonon con- 
fer the non-relativistic 



figurations. This version of the model has been proposed in Ref. 
case. It has been noticed that the energy-dependent resonant part of the two-quasiparticle 
amplitude $(a?) can be factorized to extract the two-quasiparticle intermediate propagator 
with the frequency shifted by the phonon energy. In the relativistic RQTBA it takes the 
following form: 

k 5 k 6 ,ii 

where (u — rjQ^) are the matrix elements of the two-quasiparticle propagator in the 
mean field with the frequency shifted forward or backward by the phonon energy f2 M , 

Sfc 1 fc 2 ;A;5fc6 fcifcs ' /J,;kgk2 I fj,;kik^ u kgk2 S k2ki;kaks ) V / 

and 7^ ;fclfc2 are the quasiparticle-phonon coupling vertices defined in Ref. [if]]. In the graphic 
expression of the amplitude (jSJ) in the upper line of the Fig. [T] the propagator * s 
represented by the two straight nucleonic lines between the circles denoting emission and 
absorption of the phonon by a single quasiparticle with amplitudes 7^./ Cl fc 2 - The extension of 
this model proposed in Ref. 



23| consists in the following: we introduce RQRPA correlations 
into the intermediate two-quasiparticle propagator so that the two-phonon configurations 
appear in the amplitude $(u;) as it is shown in the lower line of the Fig. HJ The analytic 



expression of the new amplitude reads: 



, j- r) j- ri* 



where 



t W =V7? (+) -V-y {+) 7? {+) =-? { - ) * (6) 

S [iv\k\k2 / j '^v\k\k§ I fJ.;k(jk2 / j I fi;kik^ '^v\k$k2 S \iv\k 2 k\i V u / 

fc 6 fc 5 



and TvJ,.*. are the matrix elements of the RQRPA transition densities defined in Ref. [If 
and corresponding to grey circles together with two nucleonic lines in Fig. [U One can 
show that in the limit of the vanishing static interaction V between the two intermediate 
quasiparticles Eq. (jSj) transforms to the Eq. ([3]) of the original (R)QTBA. As in conventional 
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(R)QTBA, the elimination of double counting effects in the phonon coupling is performed 
by the subtraction of the static contribution of the amplitude $ from the residual interaction 
in Eq. (J2J), since the parameters of the CEDF have been adjusted to experimental data for 
ground states and include therefore already essential phonon contributions to the ground 
state. Therefore, the BSE in the two-phonon model has the same form as Eq. (CQ), but 
contains the amplitude $ instead of $. In order to take the quasiparticle-phonon coupling 
into account in a consistent way, we have first calculated the vertices 7^ of this coupling 
and the transition densities W 1 within the self-consistent RQRPA using the static residual 
interaction V, as described in Ref. 15J]. On the next step the BSE for the correlated 
propagator R^ e \uj) 

R {e \uj) = R(°\lu)+R (0 \uj)1§(uj) - $(0)]I&\w), (7) 

has been solved in the Dirac-Hartree-BCS basis. Then, the BSE for the full response function 
R(u) 

R(u) = i? (e V) + R^(u)VR(u) (8) 



has been solved in the momentum- channel representations. The details are given in Ref. [15]. 
To describe the observed spectrum of a nucleus excited by a weak external field as, for 
instance, an electromagnetic field P, the microscopic strength function S(E) is computed 
as: 

S(E) = Km Q Im[Tr(P*R(E + iA)P)]. (9) 

To illustrate the effect of the additional static correlations on spectra of nuclear excita- 
tions, we consider the dipole response of tin and nickel isotopes in the area below the giant 
dipole resonance (GDR). Fig. [2] displays the dipole strength function in 116 > 120 Sn calculated 
within the conventional RQTBA [15( and the two-phonon RQTBA-2 presented here. The 
strength functions obtained in this way are compared with the original RQRPA strength 
function because both of them originate from RQRPA by different fragmentation mech- 
anisms. The first observation is that the total strength B(E1) f below the neutron 
threshold is reduced in the RQTBA-2. For example, for 116 Sn we have 0.20 e 2 fm 2 below 8 
MeV which agrees with 0.204(25) e 2 fm 2 obtained in the experiment of Ref. j^ . For 120 Sn, 
if we include the relatively strong state at 8.08 MeV into the integration region, this quantity 
is 0.31 e 2 fm 2 , in agreement with the results of the quasiparticle phonon model (QPM) of 
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FIG. 2: Low-lying dipole spectra of 116,120 Sn calculated within the RQRPA (dashed curves), 
RQTBA (blue solid curve, panels (a,c)) and RQTBA-2 (red solid curve, panels (b,d)). A finite 
smearing parameter A = 20 keV has been used in the calculations. The inserts show the zoomed 
pictures of the spectra below 8 MeV. The arrows indicate the neutron thresholds. 

0.289 e 2 fm 2 {26]. One can also notice that in both nuclei the major fraction of the RQRPA 
pygmy mode shown by the dashed curve is pushed up above the neutron threshold by the 
RQTBA-2 correlations. 

The energies and the corresponding B(E1) f values of the lowest 1~ states in the tin 
isotopes are listed in Table [0 The experimental energies, B(E1) f values and the 
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TABLE I: The energies, reduced transition probabilities and anharmonicities of the lowest l - states 
in 116 ' 120 Sn isotopes calculated in the relativistic two-phonon model, compared to data of Ref. 
see text for details. 





0,(17) 


B(E1)| 






(MeV) 


(1(T 3 e 2 fm 2 ) 




RQTBA-2 


2.66 


14.5 


0.94 


116 Sn Exp. 


3.33 


6.55(65) 


0.94 


RQTBA-2 


3.22 


16.9 


0.95 


120 Sn Exp. 


3.28 


7.60(51) 


0.92 



values, R^j = uj(1 1 )/(u(2i) + ^(3! )), are taken from Ref. [l9j. Notice, however, that the 
measurements with the larger end point energies for the electron bremsstrahlung result in 
the larger B(El) T values: 16.3(0.9) and 11.2(1.1) for 116 Sn and 120 Sn, respectively 0,0. 
A good agreement with the data is obtained in spite of the fact that these tiny structures at 
about 3 MeV originate by splitting-out from the very strong RQRPA pygmy states located at 
the neutron thresholds, that is described by the fully consistent inclusion of the two-phonon 
correlations without any adjustment procedures. 

In the RQTBA-2 the position of the first 1~ state is basically determined by the sum of the 
energies of the lowest 2 + and 3~ phonons. From the Eq. ([5]) one can see that the amplitude 
$(co>) consists of the pole terms with the poles at the energies which are sums of the two 
phonon energies. Therefore, the energy of the first 1~ state is approximately equal to the 
sum of the energies of the lowest 2 + and 3~ phonons with some small correction introduced 
by the static residual interaction V, in agreement with data 19J. This correction is given by 
the quantity R^, whose deviation from unity characterizes the two-phonon anharmonicity, 
see Table HI For example, in 120 Sn the energies of the 2f and 37 phonons calculated within 
the RQRPA are obtained at 1.48 MeV and 1.90 MeV, respectively, explaining the position 
of the I7 state at 3.22 MeV which is identified experimentally as a member of the 2f <S> 37 
quintuplet. 

The electric dipole strengths in neutron rich 68 7o,72j^j i so topes are displayed in Fig. [31 
For all three isotopes we found a redistribution of the low-lying strength as compared to 
the RQTBA calculations of Ref. jl6J. The calculated strength distribution in 68 Ni has 



8 




6 8 10 

E [MeV] 



12 



6 8 10 

E [MeV] 



12 6 8 10 12 

E [MeV] 



FIG. 3: Low-lying dipole spectra of 68 > 70 > 72 Ni calculated within the RQTBA-2 with a smearin g of 



10] 



20 keV (thin curves, pink) and 200 keV (filled area). Panel (a) contains also the data from Ref. 
(open circles, units on the right). The arrows show the neutron thresholds. 

its maximum at 10.30 MeV and the total strength below 12 MeV is 2.73 e 2 fm 2 , while 
the corresponding fraction of the energy weighted sum rule (EWSR) is 7.8% of the total 
integrated photoabsorption cross section and 11% of the Thomas- Reiche-Kuhn sum rule. 

n 

These characteristics are in agreement with the recent data of Ref. [10|]. The RQTBA-2 
dipole strength distributions presented in Fig. [3] for the 70 ' 72 Ni isotopes can be suggested as 
predictions for possible future measurements in these nuclei. 

In summary, the two-phonon version of the relativistic time blocking approximation is 
presented. Within this model it has been shown how the RQRPA modes are fragmented 
due to the coupling to two-phonon configurations, thus explaining the physical connection 
between the pygmy dipole mode and the 1~ member of the 2f <g> 3^ quintuplet. A very 
reasonable description of the lowest 1~ states in 116 > 120 Sn has been achieved within the 
fully consistent scheme avoiding any adjustment procedures. The resulting low-lying dipole 
spectra in 116 ' 120 Sn are compared with conventional RQTBA calculations. It has been found 
that the two-phonon correlations redistribute the fragmented strength as compared to the 
2q®phonon RQTBA so that the major fraction of the RQRPA pygmy mode is pushed 
above the neutron threshold and therefore mixed with the giant dipole resonance tail. The 
calculated low-energy fraction of the electric dipole strength agrees also very well with the 
available data for the tin isotopes and for the recently investigated neutron-rich nucleus 
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68 Ni. In general, the relativistic two-phonon model presented here provides a new quality 
of understanding of mode coupling mechanisms in nuclei. The method is based on Green's 
function techniques and can be widely applied also in other areas of quantum many-body 
physics. 
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